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I. INTRODUCTION 



Contemporary mathematical literature assigns an abstract meaning to the concept of 
dynamical systems. They are most often considered as substrates for particular problems 
in applied mathematics. Physics, by contrast, prefers thinking of dynamical systems as be- 
ing macroscopic objects with time-dependent properties, consisting of microscopic particles 
responding to the laws of mechanics. The goal here is to establish links between observ- 
able macroscopic properties and the more or less predictable motion or trajectories of the 
system's elements. By quoting formally fluid dynamics and statistical physics in the work- 
shop's announcement as topics to be covered, the physicist's definition is implicitly included 
in the general debate. This implies possible confrontation with the experimental reality for 
validation of the conclusions. The present review has the latter aspect of the problem of 
dynamical systems in mind. 

The general discussion concerns time-reversal properties. With physical systems, not con- 
sidering (or excluding) possible external action, the microscopic laws of mechanics governing 
the motion of the particles are Hamiltonian and therefore time-reversible. Macroscopic time- 
dependent properties show however a high degree of irreversibility. Apparent contradiction 
between the symmetry properties of the microscopic laws and the experimentally verified 
macroscopic behaviour of dynamical systems remains a controversial subject of debate. 

Macroscopic time- dependent phenomena shown by physical systems include relaxation 
of disturbances, transport effects like friction and others and bifurcations to more or less 
ordered dissipative structures. The last item has been quoted explicitly in the list of dy- 
namical phenomena to be covered by the present general discussion. Considering the many 
connections prevailing between this item and the other time-dependent effects mentioned 
above, a general discussion including relaxation and transport effects is justified. 

In this review it will be proven that next to the internal Hamiltonian laws of motion, 
observable physical systems rely on interactions with the external world for their global 
dynamics. Depending on the conditions, the step that controls the relevant time- dependent 
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phenomenon and defines the observed degree of time symmetry may depend either on inter- 
nal or external effects. In all cases a correct formalism is required for expressing the prop- 
erties of the surroundings and their relation to that of the system itself. This is achieved 
by extending the principles of thermodynamics to conditions removed from equilibrium. A 
review of paradigms and assumptions supporting the time-reversal symmetry problem in 
statistical physics will be given first. They will be critically analyzed and a new approach 
will be suggested based on analysis of Joule's experiment. The necessary thermodynamic 
tools will then be developed. The thermodynamic formalism will be applied for predicting 
and discussing selected topics in fluid dynamics. In a first step, the transport coefficients 
are considered. Structure formation and transition to turbulence are developed from there 
on. 



II. CONTROVERSIES 

Distortion of macroscopic systems initiates spontaneous and irreversible processes tend- 
ing to restore the previous state of equilibrium or possibly to establish a new one. This fact 
is most often taken to be a genuine property of conservative Hamiltonian systems. Conflict 
between the time asymmetric behaviour of macroscopic relaxation and the strict reversible 
nature of microscopic dynamics of conservative systems has been the subject of discussions 
over nearly a century [§]|§. 

Results published in recent decennia in applied mathematics concerning time dependent 
transformation of systems where the number of identical elements tends to infinity at con- 
stant density HU have been a stimulus for trying to solve the irreversibility paradox. The 
arguments are related to the mathematical property of mixing. This is meant to express dis- 
semination of the system's parameters throughout the available space towards homogeneous 
and statistically independent distributions. It associates irreversibility to an infinite Poincare 
recurrence time for most initial fluctuations. Furthermore, progress obtained in characteriz- 
ing deterministic chaotic motion (Lyapounov exponents) spurs theoretical research towards 
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relating the relevant numbers to transport properties associated with irreversible relaxation 
dynamics. In this context chaotic scattering of particles on hard disks (Lorentz gas) is 
used as a model for describing irreversible diffusion in agreement with reversible microscopic 
dynamics 

Using the language of functional analysis one arrives at a rigourous formulation of the 
fundamental difference between reversible and irreversible evolutions. The reversible evo- 
lution of dynamical systems is expressed by a group of isometric operators, whereas the 
irreversible one is described by a semigroup of strictly contractive ones. As microscopic laws 
are time-reversible, the evolution they cause must be formulated by isometric groups. By 
contrast, macroscopic irreversible processes can be described only by contractive semigroups. 
We are thus faced with the question of how this transition from a group to a semigroup can 
be justified. The presently most popular answer is proposed by the Brussels school associ- 
ated with the name of I. Prigogine ||. This school's interpretation rejects formally the role 
of loss of information in deriving the second law of thermodynamics. It has been subjected 
to criticism in [|T(J] on the basis of mathematical results of G. Braunss [[TT[ . 



Aforementioned research trends assume explicitly that irreversible processes would occur 
in systems isolated from the environment. The statistical properties of the time dependent 
random forces acting on the particles are supposed to be completely determined by the initial 
conditions and by the dynamics of the system ||. Rescaling is said to make phenomenological 
parameters (e.g. viscous drag controlling the motion of Brownian particles) converging 
to genuine irreversible properties of the system. Dissipative coupling with a reservoir is 
therefore explicitly rejected, as "artificial and unnecessary" @||. 

Contrasting with the latter, some authors still insist on the unavoidable interaction of 



macroscopic systems with their environment acting as a reservoir or heat bath ||T2| , |T3|| . Stress 
is laid on the environment where the fluctuating forces to be introduced into the equations 
of motion of the system of interest come from. Dissipation arises from back-reaction of the 
environment to the evolution of the system [0,0 • As an example, concerning viscous drag 
as mentioned above, the relevant dissipative force is said to originate from interaction of the 



Brownian particle with the surrounding fluid acting as a reservoir. The time averaged value 
of this force gives rise to Stokes' law, on which fluctuations are superimposed. 

To anticipate and prevent future ambiguities, we specify the meaning to be given to some 
important key-words frequently used in the discussion to follow. 

By isolated systems are meant systems the dynamics of which is defined in a unique 
way by time-independent Hamiltonians. They are obviously conservative. Since we are 
focussing next on systems that are not perfectly isolated, we need to give them an appropriate 
name. A system will be called quasi-isolated if it exchanges energy but not matter with its 
surroundings. 

Equilibrium refers to any macrostate that is stationary (i.e. time-independent) under 
the given macroscopic constraints. As an example we may consider that this property holds 
for a gas contained within walls at different temperatures when it has reached stationary 
conditions. By approach to equilibrium we denote any observable macroscopic evolution 
towards equilibrium. 

III. ISOLATED DYNAMICAL SYSTEMS 
A. The Problem 

Consider a system of N identical structureless classical particles constrained to move in 
volume V. The state of the i-th particle is in unique way given by the vector of generalized 
coordinate q\ and by the vector of generalized momentum pi. Expression (q,p) will denote 
the state of the collection of iV particles, or a point in the system's phase space T. This 
system's dynamics is described by Hamilton's equations with the time-independent Hamil- 
tonian H(q,p) ||16| . The action of the walls is represented by a high repulsive potential (a 
possible choice being an infinite potential barrier). The solution (q(t),p(t)) to these equa- 
tions for some initial state (g(0),p(0)) can be written using a group of dynamical evolution 
operators T t : 
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(q(t),p(t))=T t (q(0),p(0)). 



(1) 



It represents the motion of a point along a trajectory in phase space T. 

Investigation of the evolution of an individual system is apparently not possible and 
therefore one resorts to statistical methods fT7| . One considers an ensemble of equivalent 
systems. The ensemble properties are then described by means of a distribution function 
f(q,p] t), giving the probability density of finding a system at point (q,p) in phase space at 
time t gj|. 

From the deterministic description defined by the group of operators T t , one goes over 
to the statistical description of evolution of the distribution function by introducing a group 
of operators U t defined through relation 



It is easy to show that the statistic evolution operator Ut is unitary — and thus isometric 
- if it acts in a Hilbert space C 2 (T) of square- integrable functions defined on phase space 
T Jil|]. By Stone's theorem there exists then a self-adjoint operator L (Liouville operator), 
such that U t = e~ lLt . Using the Poisson bracket notation it can be easily proved that L can 
be expressed by means of the Hamiltonian function as Lf = i{H, /}. 
The Liouville equation, 



has for solution: f(q,p;t) = U t f{q,p;0) = e tLt f{q,p;0), where f{q,p;0) = /(g(0),p(0)). 



Let if denote one of the possible equilibrium functions. Since if does not change with 
time, dip/dt = 0, whence Lip = 0. From (^) we have obviously U t f = f- Now, unitarity of 
U t (i. e. the equality Z7 t * = U-t) implies, first of all, invariance of the scalar product in the 
course of evolution 



This reduces to = ||g|| = const for / = g. Taking now g = f — f, where / is any 

nonequilibrium distribution function, we obtain 



Utf(q,p) = f(T_ t (q,p)) . 



(2) 




(3) 



{U t f | U t g) = (f\g) for any f,g£ L 2 (T) . 



(4) 
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Utf — ip\\ = const . 



(5) 



Hence, approach of mechanical systems to equilibrium cannot be described statistically 
as a strong convergence of the distribution function to an equilibrium function. 



The description of a system on the basis of a distribution function f(q,p;t) is complete 
from a statistical point of view. However, solving equation (|3|) explicitly is apparently 
not possible. Since the distribution function cannot converge strongly to any equilibrium 
function, the entropy as defined by Gibbs by relation 



is equally time- independent ||17| . This drawback is eliminated by introducing so-called re- 
duced distribution functions which are obtained by integrating function f(q,p;t) over a 
subset of particles. For example, we obtain an r-particle distribution function f r (q,p;t) by 
integrating f(q,p; t) over positions and momenta of iV — r particles (0 < r < N — 1). Such 
a function describes statistical properties of a subsystem comprised of r particles. 

It turns out that subsequent integration as above of Liouville's equation (fj) leads to a 
system of N—l integro-differential equations for a set of N— 1 reduced distribution functions. 
The equations have the property that the equation for f r (r < N) contains always function 
f r+ i too |p0[ . This system (or hierarchy) of equations, the BBGKY hierarchy, can be solved 
only if interrupted at some level. To break the hierarchy, say at the r- th level, means 
expressing function / r+1 with the help of f r in equation for f r . Such an interruption having 
been performed, the equation for f r is in principle solvable. The solution may subsequently 
be substituted in the equation for / r _i. The procedure may be continued until we are left 
over with the only equation for the single- particle distribution function fx. This is the 
master (kinetic) equation. 



B. BBGKY Hierarchy 




(6) 
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The mathematics leading to an f r independent of term f r +i, can obviously not be sub- 
stantiated rigourously. One relies exclusively on more or less acceptable heuristic semi- 
quantitative physical assumptions. 

The first attempt in this direction was proposed in 1872 by L. Boltzmann, who derived 
an integro-differential equation for the one-particle distribution function of a gas and proved 
that its solution proceeds to the equilibrium distribution function as time goes on. Among 
the assumptions used, the most important one was the so-called molecular chaos hypothesis, 



according to which particle momenta are independent of their positions This hypothesis 
is eventually equivalent to the assertion that the evolution of a system is not fully deter- 
ministic and may be regarded as a stochastic process. Hence the derivation of second law 
of thermodynamics (or H-theorem) from Boltzmann's equation cannot be considered as the 
proof for approach of gases to equilibrium. Despite this, Boltzmann's equation has become 
the foundation of the kinetic theory of gases. The hypothesis of molecular chaos remained, 
however, for more than a century the target of either criticism or efforts to find its exact 
justification . 



C. Spectral Theory and Approach to Equilibrium 

A well-know theorem of functional analysis [p3fl , states that the spectrum a(L) of any 
self-adjoint operator L, acting in a space £ 2 (r) of square- integrable functions, can be decom- 
posed in a unique way into three constituents: the pure point spectrum o~ pp (L) consisting of 
eigenvalues of L, the absolutely continuous a ac {L) and the singular spectrum a S i ng (L). Thus 

a(L) = o-pp(L) U a ac (L) U a sing (L) . (7) 

Let {E\} be the resolution of identity for operator L and Vf = (E\f | /) the spectral 
measure generated by a function / G C 2 . Decomposition (|7p of the spectrum implies de- 
composition of C 2 in orthogonal subspaces 

= £pp ® ^ac ® ^sing > (8) 
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where C pp = {/ G C 2 : Uf is pure point}, £„ c = {/ G C 2 : z// is absolutely continuous} and 
C 2 ing = {f E C 2 : Uf is singular} p3| . The basis of C 2 p consists of eigenf unctions of L. We 
shall denote / G C 2 pp by f pp , f G C 2 ac by / ac and / G £^ ng by / sing . Decomposition (§) means 
that every / G C 2 can be written in the form 

/ fpp ~l~ fac 4" fsing ■ (9) 

It must be stressed that this decomposition is invariant under the action of evolution operator 
Ut, i. e. U t C 2 c = C 2 ac1 U t C pp = C pp and U t C 2 ing = C 2 ing . This means that each component 
of / evolves under U t independently of the other ones. Hence, the components evolve in 
orthogonal subspaces of C 2 . The dynamics of the system can be, therefore, decomposed into 
subdynamics describing qualitatively different types of evolutions. 

In view of the above, any component of a distribution function evolves with time in- 
dependently of the other ones. One can, therefore, investigate the time evolution of each 
component separately. It turns out that for realistic Hamiltonians (and therefore Liouvil- 
lians) a(L) sing = 0. Actually, the singular part of the spectrum corresponds to "patholog- 



ical" Hamiltonians P^j . For f ac it was shown in |25| that Utf ac tends weakly to zero for 
\t\ — > oo, i. e. in this limit we have (Utfac I <?) — for any g G C 2 . In an alternative notation, 
Utfac —> 0. Finally, the evolution of U t f pp has been proved to be quasi-periodic. 

Functions with the property U t f = f (for any t) span obviously a linear subspace. Let us 
denote the latter by C 2 const . For each / G £l onst we get from @ Lf = 0, i. e. / G KerL. The 
inverse implication is equally true, so that we have C 2 const = KerL. A solution of equation 



U t f = / is an eigenfunction of operator L corresponding to the eigenvalue A = 0. This 
' 2 C r 2 

'const — *"pp ' 



yields KerL = C 2 const C C 2 



Let now M C T, y,(M) = e, < e < oo and let x(M) G C 2 (T) be the characteristic 
function of the set M (i. e. x{ x ) = 1 f° r x E M and xi x ) = otherwise). For arbitrarily 
small e > 0, the convergence Utfac — * implies 

lim (Ufac I x(M)) = , (10) 

\t\^oo 



which is equivalent to 

lim / (UJ ac )dfi = . (11) 

This resembles the behaviour of function sin(tx) on the finite interval of R and one can 
regard convergence ( |TT| ) as the generalization of the Riemann-Lebesgue lemma to functions 

The above implies that U t (f ac + f const) tends weakly to f const- Weak convergence is math- 
ematically a well defined property with physically relevant contents. In physical systems, 
stress is rather lead on the observables, or expressions like (Utf\g), where g represents 
dynamical or generating functions. Hence, if the Liouville operator acting in C 2 (T) has 
an absolutely continuous spectrum and if its only eigenvalue is A = 0, then, for positively 
definite functions (distribution functions), as \t\ — > oo, (U t f \ g) tends to (U t f \ f CO nst), that 
is the value of the observable corresponding to an equilibrium state. This result may be 
regarded as the exact criterion of approach to equilibrium in conservative systems. 

The above results show that approach to equilibrium in isolated systems makes sense 
only for observables, not for distribution functions. Distribution functions themselves change 
reversibly, and without any apparent tendency to approach equilibrium. Observables related 
to different distribution functions may tend to the same limit, meaning that the approach to 
equilibrium is connected to a loss of information caused by the integration leading to (Utf \ g). 

D. The Brussels School Theory 

As mentioned above, the spectral theory of operators shows that the spectrum of any 
self-adjoint operator L is real and that operators defined as U t = e~ lLt , with self-adjoint L, 
represent a one-parameter group of unitary operators preserving the norm, i. e. ||f/t/|| = ||/|| 
for any / e C 2 . Operators preserving the norm are called isometric. 

Unitary group Ut is defined for any real t and — if such a group describes the evolution of 
distribution functions — this means, if an evolution is possible for t > 0, then the evolution 
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for t < is equally possible. Hence, the unitary group describes a reversible evolution, 
expressed explicitly by U_ t {U t f) = f. 

Irreversible evolution can be described only by semigroups (which cannot be continued 
into groups), i. e. sets of operators defined e. g. only for t > 0, not for t < 0. Let us 
denote such a semigroup by W t . The impossibility to define a semigroup for any real t 
follows from the spectral properties of its generator K: Wt = e~ lKt . If the spectrum of K 
contains points with negative imaginary components, the semigroup is strictly contractive, 
i. e. ||Wi/|| < ||/||, so that it is not isometric and it cannot be continued into a group |26|| . 



The reversible evolution is expressed by a group of isometric operators, the irreversible 
one, on the contrary, by a semigroup of strictly contractive operators. The most popular 
answer to the question of how such a transition from a group to semigroup occurs when 
going from the microscopic level to the macroscopic one is given by the Brussels school. The 
short formulation of this viewpoint might run as follows. Suppose that for a given system 
there exists a nonunitary operator A such that 

(a) W t = MJ t hr l is a strictly contractive semigroup for t > 0, 

(b) W t is positivity preserving, i. e. f{uj) > for almost all u G T implies W t f(uj) > for 
almost all u too, and 

(c) WtJ "const f const' 

The dynamics U t of the system is then said to be inherently stochastic, meaning roughly 
that the microscopic behaviour of the system is deterministic, while it behaves as if it was 
really stochastic. That is why it is appropriate to describe the dynamics using statistical 
methods [27|. Under the above conditions it is claimed specifically 0, that operator A 



converts the deterministic evolution into a Markov process whose transition probabilities 
can be obtained from the adjoint of W t by 

P(u,A;t) = W; X (A) , (12) 

where P(uj, A; t) is the probability of transition from point to to domain A of phase space in 
time t and x(A) is a characteristic function of the set A. Operator A is called the system's 
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Lyapunov converter. 



It is not expected that A would exist for all dynamical systems (flows). It is shown in |28 
that the mixing property is necessary and the condition of if -flow sufficient for the existence 
of a Lyapunov converter A. The close connection between the intrinsic irreversibility of the 
system, expressed by the existence of A, and the instability of the motion is put forward. 

The result of the Brussels school is important in that it divides the systems into two 
categories: systems in which it makes sense to expect thermodynamic behaviour and systems 
where this is not the case (trivial examples of the latter are harmonic oscillators). This 
classification is however used as the basis for derivation of far-reaching consequences which 
are not always convincing. 

The Brussels school does not give any physical interpretation to the action of operator 
A. It is only stated that stochastic evolution arises from a deterministic one "simply as a 
result of 'change of representation' brought about by (non-unitary) similarity transformation 
A" |J. The physical meaning of this transformation is indeed nowhere explained. Moreover, 
the analysis of the properties of Markov processes shows that the claim, that every semigroup 
W t of operators acting on C 2 and having the properties (a)-(c) (see above) defines a Markov 
process, is not in general sufficient to prove that the process has time-independent transition 
probabilities. This condition is necessary for a process to have any physical meaning at all. 

The action of A is equivalent to a change of dynamics of the system from deterministic 
to stochastic [p9|j . Such a change can be performed in many ways and the operator A is 
constructed precisely in such a way that the transition probabilities do not depend on time. 
The possibility to change deterministic evolution into a stochastic one does not give any 
information about the original deterministic evolution, as is demonstrated by the following 
example fTO 



Consider the unitary group of shifts U t f(x) = f(x+t), f G C 2 (R; dx). Using the operator 
of differentiation d/dx, this group may be written, in the form 

U t = e td/dx , (13) 
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so that the generator of the group is operator L = id/dx. For this group it is possible to 
construct a Lyapunov converter A, defined on the subspace of continuous functions with 
compact support, converting U t into the strictly contractive semigroup W t = exp(t d 2 /dx 2 ), 
t > associated with a one-dimensional diffusion process, like diffusion of heat [[□]] . The 
group of shifts U t f(x) = f(x + t) lacks obviously any resemblance to a process that could 
be either random, stochastic, approaching equilibrium, or even irreversible. The only con- 
vergence property it has is the weak approach to a function ip = 0. 

This example shows clearly that the existence of a Lyapunov converter (be it that con- 
verting an isometric group into a strictly contractive one) is not sufficient to prove the ap- 
proach to equilibrium. Consequently (in view of our comments about the relation between 
irreversibility and approach to equilibrium) this cannot justify the origin of irreversibility. 
However, it turns out that the notion of Lyapunov converter is significant in understand- 
ing the essence of the method of complex scaling introduced originally for identifying the 
so-called resonant discrete points in the continuous spectrum of quantum systems pOj-pTJ . 



IV. QUASI-ISOLATED DYNAMICAL SYSTEMS 

The discussion above was centred on the quest for solutions to irreversible approach to 
equilibrium on the basis of reversible microscopic laws in perfectly isolated systems. The 
proposed theories leave us disappointed, being restricted to qualitative results and concerned 
only with the limiting behaviour for t — > oo. They are unable to predict finite-time dynamics 
and experimentally verifiable results. 

Contrasting with the above let us not consider the system's boundaries (walls) as a 
potential to be added to the Hamiltonian, but instead as the locus for an independent 
contribution to the particles' motion. This is the quasi-isolated system's paradigm. 
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A. Experimental evidence 



Discussion of the directionality of time's arrow is often introduced intuitively on the basis 
of a simplified representation of Gay-Lussac's experiment. A box is considered, consisting 
of two compartments, the parts being filled with gas at different pressures. Prior to the 
experiment the gas is assumed to be at equilibrium. The long time evolution towards a new 
equilibrium distribution following rupture of the division is taken to be modelling irreversible 
behaviour of the global dynamics. Joule repeated Gay-Lussac's experiment with great ac- 
curacy. His purpose was to measure possible heat exchange with an external calorimeter 
associated with spontaneous expansion. With an ideal gas, if no mechanical work is allowed 
to be performed during the process, when the system had reached its final state of equilib- 
rium, no net exchange of heat with the surroundings was observed. Joule concluded that 
the system behaved as if it were isolated. Be it stressed that Joule was considering only the 
initial and the final conditions, neglecting whatever dynamics was involved in reaching final 
equilibrium. 

Let us make the experiment more realistic by examining the effect of pricking an air- 
inflated balloon inside either an acoustic reverberation hall or an anechoic chamber. In 
both cases the excess air contained in the balloon disseminates spontaneously throughout 
the rooms, never to come back again, compressed in its initial volume, but the subsequent 
process is very different indeed. In the reverberation hall an acoustic perturbation is created 
and, the better the walls' reflecting quality, the longer it remains. By contrast, in the 
anechoic room, the perturbation vanishes promptly. In the reverberation room, some energy 
is stored in a coherent or collective motion (acoustic perturbation) where it remains as the 
memory of the initial conditions. With correctly shaped walls, the initial information might 
even be partially retrieved as echoes. By contrast, in the anechoic room, memory of the past 
is soon forgotten. 

Initial and final conditions are identical in the two cases and so is the air inside the rooms, 
and therefore the frequency and the quality of the inter-particle collisions (Hamiltonian 
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dynamics) assumed usually to be the source of relaxation. The only difference between the 
two experiments is the nature of the walls. One is therefore forced to conclude that global 
relaxation dynamics of a spontaneously expanding gas leading to final equilibrium depends 
on the acoustic (physical) quality of the walls representing the system's environment. 

The experiment suggests convincingly that the global dynamics consists of two indepen- 
dent major steps. According to the coupling efficiency of the system to its surroundings 
(impedance matching), either step may be rate determining. If coupling is very effective, 
global dynamics is controlled by slow transport of mechanical properties to the walls (mo- 
mentum, energy, matter). Transport coefficients like viscosity, thermal conduction etc., 
are correctly defined only in such non isolated conditions. By contrast, if the system is 
nearly isolated (quasi-isolation) , memory of the initial conditions remains for some time as 
a collective or coherent motion of the particles (acoustic motion) and full thermodynamic 
equilibrium is slow to reach. Strict isolation and transport effects are incompatible. 

The two steps involved by the scenario are very different in their dynamics. Depending 
on the system of interest, they may be more or less concomitant. For simplicity, we shall 
take them next as frankly separated in the time. 

In Joule's experiment, as soon as the membrane has been ruptured, a stream of gas is 
ejected from the compartment at the highest pressure, thereby creating a collective motion 
of the particles. By performing work on itself, the system transfers energy into the jet. This 
is subtracted from the initial thermal supply (adiabatic expansion). Loss of thermal energy 
is equivalent to cooling. 

On reaching the wall opposite the puncture, if this is hard, the initial jet is reflected 
and turns progressively into a compound acoustic perturbation with the same energy. The 
spectrum and phases of its components are the memory of the initial conditions and of the 
shape of the reverberating walls (coherence). This is the first step of the general process. 
When this is done, although the particles are disseminated throughout the whole volume, 
the system may not be claimed to be at equilibrium (weak irreversibility) . 

Relaxation of the coherent or collective motion starts now. Every collision with inco- 
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herently fluctuating wall atoms interrupts the running canonical trajectory. A new one 
starts, with possibly modified initial conditions^. Earlier correlations are progressively bro- 
ken, thereby thermalizing the energy accumulated initially in the jet and stored later in 
the acoustic perturbation. As a result, if the system is an ideal gas (hard spheres allowed), 
aforementioned transient cooling is progressively neutralized, as is expected by Joule's re- 
sult. When final equilibrium has been reached, collective or coherent motion has relaxed 
and information about the initial conditions is completely lost (strong irreversibility) . 

Initial dissemination of the particles throughout the system follows conservative Hamil- 
tonian dynamics. No matter how intricate (chaotic) the motion of individual particles may 
be (Sinai billiards), this part of the motion preserves the memory of the initial conditions. 
Contrasting with the latter, the correlations removing step, where the particular properties 
of the walls determine how efficiently the system is coupled to its surroundings introduces 
in the global motion stochastic non Hamiltonian jumps between the different accessible tra- 
jectories. According to whether during the particular impacts the relevant wall atom moves 
towards the colliding particle or in the opposite direction, transient work is transferred to 
the system or to the environment. Energy fluctuates about its average value. Only if the 
thermodynamic requirement is fulfilled that the system and its neighbourhood are at the 
same temperature does the average energy transfer vanish. 

B. Thermodynamics 



x It should be obvious that this behaviour is much less stochastic than the one put forward by 
Boltzmann's Stofizahlansatz, according to which stochasticity builds up at every collision between 
two particles. 
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1. Entropy 



Any function determined completely by the set of constraints defining the relevant sys- 
tem's particular macrostate is a function of state. In 1865 Clausius discovered a function of 
state called entropy which, for reversible processes, is defined as a differential SS = SQ/T, 
where T is the system's temperature. In 1877, Boltzmann derived an expression that links 
the experimental entropy to the statistical properties of the relevant macroscopic system. 
This reads 

S = k B ln[W(A)], (14) 

with meaning the probability associated to macrostate A. The latter is to be inter- 

preted as the total volume accessible to the motion in phase space, given the set of constraints 
(represented by the collective variable A) describing the system's particular macrostate. Let 
it be noted that equilibrium macrostates are usually defined by their total energy E, particle 
number of any sort N r and physical volume V, the traditional microcanonical variables. In 
the literature, extension of the discussion to non-equilibrium macrostates is avoided. This 
limitation will be reconsidered below. 

In a strictly conservative isolated environment, the dynamics being described by a sin- 
gle multi-particle trajectory in phase space, no matter how intricate (chaotic) this may be, 
transitions between different trajectories are not possible. Then, according to the definition, 
the entropy is zero and it does never change. This conclusion is consistent with Liouville's 
theorem claiming conservation of the measure in phase space when the mechanics is conser- 
vative. 

For Boltzmann's entropy to be a pertinent function of state, prompt accessibility of all 
the quantum states or trajectories belonging to the given macrostate is required. Acces- 
sibility means incoherent transitions between the available and accessible trajectories or 
quantum states during the observation period. This depends on fast uncorrelated action of 
the environment with fluctuating exchange of mechanical properties (momentum, energy). 
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As a corollary, and as expected by the statistical nature of the thermodynamic functions, 
it appears that the definition of the entropy implies some averaging over the time. The 
resolution linked to the definition of the entropy is the average lifetime of conservative tra- 
jectories. With macroscopic systems, where the impact rate with boundaries goes to infinity, 
the average lifetime and hence the time resolution tend to zero. 

Relaxation implies relief of constraints. It opens the way to an enhanced choice of 
trajectories (microstates). Accessibility of more trajectories increases Boltzmann's entropy. 

In describing equilibrium states, the extensive variables mentioned traditionally are the 
basic microcanonical constraints E, V and N r . In order to specify unambiguously non- 
equilibrium macrostates, where more constraints prevail, additional extensive properties 
must be included. This may be for example the momentum associated with a possible 
collective or coherent motion of the system, where some of the total energy is stored (e.g. 
the jet or the acoustic motion in the aforementioned Gay-Lussac experiment). Many other 
possible distortions with respect to equilibrium may occur, like moments of the energy or 
density distribution, etc. 

Let the list of the extensive properties defining the constraints of a macroscopic system 
in a particular macrostate be written {Xi}. The entropy is a function of this collection of 
variables. By differentiating the entropy with respect to the set we get by definition the set 
of conjugate intensive variables or intensities {£/}: 

dS = E 7TF dXl = ~ kB E & dX i- ( 15 ) 
i dA i i 

This equation may be considered as defining the temperature (dS/dE) -1 and the chemical 
potential —T(dS/dN r ). In non-equilibrium conditions it generalizes all the definitions by 
proposing an intensity conjugate to each of the additional non-equilibrium constraints. 

Equation fll5l) is Gibbs' celebrated equation, generalized to non-equilibrium macrostates. 
In the simplified model of a spontaneously expanding jet mentioned above (velocity of the 
collective motion is v), the new version of Gibbs' equation reads 

dS = ^ + ^dV-Y: |diV r -k B a- dP, (16) 
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where P = Nmv represents the collective momentum of the jet and a the conjugate intensity. 
It may be shown |15[ that a = v/ksT. In the last term of equation flilf ), the differential of 
the collective or coherent energy is easily recognized. We have therefore equivalently 

AS = + ^dV — ^2 TjrdN r — — d (coherent energy). (17) 

Energy conservation throughout the expansion makes dE = and isolation causes 
dN r = 0. During the adiabatic dissemination period, the second term (work made available 
by expansion) is very exactly balanced by the last contribution (energy stored in the coher- 
ent motion), making dS = 0, in agreement with Liouville's theorem for isolated conservative 
motions. Final relaxation involves transformation of the coherent motion into thermal en- 
ergy. When this has been achieved, thanks to stochastic exchange at every impact with the 
boundaries, the integral of the last term vanishes and Gibbs' equation yields the correct 
final equilibrium entropy after expansion. 

If the contribution regarding the non-equilibrium constraint had been omitted in equa- 
tions flT6| ) and (|TTD, we would not have been able to describe the thermodynamics of the low 
entropy non-equilibrium transient state. 



2. The generalized Massieu function 

The inconvenience of considering the entropy as the leading thermodynamic function is 
that it is an explicit function of the extensive properties (Xi), while the intensities are 
better measured and controlled by the environment. That is why thermodynamics makes 
widely use of potentials and other Massieu-Planck functions, obtained from the entropy or 
the energy by Legendre transformations. 

Most popular are free energy transformations. However, contrasting with the second 
law concerning the entropy, general laws involving the energy do not exist. It is therefore 
advisable to consider transformations involving the entropy itself. If all the parameters 
(excepting the system's physical volume V) defining non-equilibrium conditions are included 
in the transformation, we obtain the generalized Massieu function V): 
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M(&,V) = £ + £GX,. (18) 
Kb , 

Unlike Massieu's original proposal, .M is an explicit function of all the state defining inten- 
sities. It may be verified that 

The advantage of referring to a state function depending explicitly on intensities is that, 
with promptly exchangeable properties, the relevant intensities of the system of interest re- 
main at all times equal to their values in the neighbourhood. We might call them strong in- 
tensities (e.g. the temperature in an efficiently thermostated system). Dynamics of transient 
states refers to the intensities of the rate determining slowly exchanging or soft properties. 

The Legendre transformation changes the maximum entropy condition with respect to 
fluctuations of the extensive variables into a minimum of the generalized Massieu function 
with respect to the intensities relating to non-exchangeable properties (e.g. particle numbers 
and their distribution in closed systems). With transient effects, this fundamental property 
defines the path followed by the system during relaxation. It gives a key for treating coupled 
flows. 

The equations above (|I8"|- [I9D are generally valid. With ideal gases, the expression for the 
generalized Massieu function takes a very simple form. Individual motions being indepen- 
dent, the global motion may be represented by a swarm of points in a reduced 6-dimensional 
single-particle phase space From here on, T will represent a one-particle phase space. 

Let f(x), x G T, be the most probable particle distribution: that which maximizes the 
entropy. Any extensive property Xj is then related to a generating function 4>j(x) so that 



x 3 = Mx)f(x)ar = (hif) . (20) 



In that context, function f(x) is readily known to be (15 



/(x)=exp[]T^(x)], (21) 
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where the intensities are the Lagrange multipliers used in the maximizing process for 
the entropy. 

With the latter distribution function, it may be verified that M. takes the very simple 
form 

Mfo,V) = I expE6MaO]dT. (22) 
Jr i 

Its numerical value is the (average) number of particles contained in the system. Through 
the integration limits in configuration space it has the system's physical dimensions (V) 
as one of its independent variables. By restricting the integration to the only momentum 
coordinates, a local generalized Massieu function is obtained, the value of which represents 
the average local density in configuration space. 

With real gases, the generalized Massieu function is modified due to the interaction 
potential between the particles. The simplified formulation is however still useful as a 
low density approximation when the duration of the inter-particle collisions is negligible 
compared to the time separating the collisions (e.g. hard spheres). 



V. TRANSPORT COEFFICIENTS 

The main objective of the theory is to predict transport coefficients from first principles, 
to be compared with the experiment. Since Boltzmann's equation there has been a consid- 



erable literature concerning that question p0| . |32| . Most frequently cited are the traditional 



Chapman and Enskog derivations |33| and the more recent Green-Kubo formalism [^(J . 

When referring relaxing systems to their fixed boundaries, the thermodynamic approach 
fits best into an Eulerian frame [j3~4 , [35 : 



^ = {f N ,H} + J. (23) 

Poisson bracket {/W, H} expresses implicit deterministic contribution to the motion while the 
source/sink term J describes explicit stochastic action of the environment. If df^/dt ^ 0, we 
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have a transient state. Either implicit or explicit contributions may then be rate determin- 
ing. If d/jv/dt = 0, J may still be different from zero, balancing a non-zero {f^,H}. That 
are stationary states. Steady transport of heat between reservoirs at different temperatures 
and steady transport of momentum in the Couette flow belong to that class of processes. 
Stochastic action of the neighbourhood defines a finite life-time to deterministic trajectories. 
In stationary conditions, flows of extensive properties supported by J are obtained by inte- 
grating the deterministic contribution over this average life-time. For individual particles, 
this is the average periodicity r of effective relaxing collisions. 

In very low density systems, when the mean free path is comparable to the system's 
physical dimensions (Knudsen gas), properties picked up by the particle at one wall are 
transported in a single jump to the opposite one. In the thermodynamic limit (non-Knudsen 
regime), head-on collisions of identical particles do not hamper transport properties but 
parallactic or off- axis inter-particle collisions do. They reduce the range of free transport, 
while information about the reservoir conditions is transferred to the relevant bulk region. 

We consider a given extensive property Xj with generating function <f>j(x), x G T, and 
we investigate its flow along the z- direction. Let z* be the ordinate of an arbitrary plane. 
The basic equation for the relevant flow Jj through this plane is 

J i = \ I I I ^ f , , ^ «pE&M*)1 dz. (24) 

T J J J a J(z*-p z T/m) £ 

The term 1/h 3 is introduced for the sake of normalization [ j36fl . 

Equation (|24]) stresses that free transport of the given property by particles is limited 
to the life-time of their trajectories. Local thermodynamic conditions at the latter's onset 
determine how much of the property is transported. Hence, the lesser the collision frequency, 
the more effective is the transport. Collisions increase the system's resistance to flow. 

When integrated, the effective collision frequency comes in the normalized form 
r/{Dy/f3rn). With Knudsen systems (mean free path at least of the order of the system's 
physical dimensions), this parameter equals 1. In thermodynamic conditions, the parameter 
is much less than 1, justifying expanding the integrand to its lowest order in r. 
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In the following, equation ( |24| ) will be applied to different types of flows. 

A. Single-component gases 

1. Viscosity 

We consider a fluid bound by a pair of walls distant by 2D, moving in opposite directions 
(Couette flow). The system's stationary conditions are defined completely by the set of 
constraints listed in table |. The intensities under direct control of the surroundings are 
the particles number, the kinetic energy and the intensity conjugate to the gradient of shear 
momentum. It may be verified [jnj that the velocity of the walls (y— direction) equals ±a y /f3. 



TABLE I. Constraints for Couette flow (( = z/D) 

* Mn 

Particles number 1 a 

2nd moment of particle distribution [£ 2 — 1] 62 

Kinetic energy ^(p 2 /2m) — (3 

2nd moment of energy distribution (£ 2 — 1) ^(p 2 /2m) —72 

Gradient of shear momentum (p y a y 
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Two variables remain to be determined, namely 82 and 72, requiring two independent 
equations. In stationary conditions, pressure gradients or acoustic perturbation are absent. 
This is expressed by vanishing z*-dependence of flow of transverse momentum (p z )- Likewise, 
the total flow of energy through the system is zero. By implementing equation (|24] ) with the 
two relevant generating functions, the conditions dJ p Jdz = and Je = yield together 

2 = 0, ~ ? = (25) 

' 2 /3 2(3 K ' 

Flow of the shear component of momentum (p y ) is obtained by implementing equa- 
tion (|24]) with generating function p y , where 9 2 and 72 are replaced by their values. This 
yields 

j (26) 
Jp »~ 2(3D /?' 1 DJ 

where n = M./V is the particle density. 

The coefficient of shear viscosity is the ratio of the forces applied to the plates, com- 
pensating for transfer of momentum from wall to wall, to the velocity gradient (a y /(j3D)). 



Following equation (p6|) , its value is 



r) = n T -. (27) 



2. Thermal conduction 

Let us consider now a system in thermal contact with two planar heat reservoirs at 
different temperatures separated by 2D. The system's stationary conditions are completely 
described by the set of constraints listed in table |T[ By inspecting the generating function 
conjugate to the temperature gradient it is clear that fc^VT = —^i/(f3 2 D). 

In the presence of a temperature gradient, particles moving towards the cold wall have 
been equilibrated with the system upstream in a hotter region at the instant of their last 
collision and vice-versa. Hence, in moving from the hot wall to the cold one, particles travel 
on the average faster than in their return cycle. If the particles are to change their average 
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kinetic energy in a correlated fashion on impact with either walls, while the container (the 
pair of walls) is to remain on the average immobile, collective momentum must be transferred 
by the container into the system. That is why a generating function for collective motion 
of the particles perpendicularly to the walls needs to be considered in constructing the 
expression for flow of heat. 

TABLE II. Constraints for thermal conduction (£ = z/D) 



Xl 6_ 

Particles number 1 a 

Gradient of particle distribution £ 6\ 

Kinetic energy J2(P 2 /^ m ) —(3 

Gradient of energy distribution Q/D J2(p 2 /2m) —71 

Collective momentum p z a z 



25 



The intensities under direct control of the surroundings are the particles number, the 
average kinetic energy and the temperature gradient. Two intensities remain to be deter- 
mined: 9\ and a z . This requires two independent equations. One is stationarity. The other 
equation describes mechanical equilibrium of the system between its walls. 

Let z* be an arbitrary position between the boundaries. The average local density n 
responds to the equation 

n ( z *) = J^j J j ' dpxdpydp* exp{£) & <M(z = z *),Px,P y ,Pz}}- (28) 

We call n + (z*) the partial density of the particles with positive velocity along the 
z— direction. Stationarity implies that this partial density equals the sum of the densi- 
ties of the particles present at places from where they will be reaching this position without 
disturbance after one collision period, their velocities being oppositely oriented. Hence 

[ />oo />oo rO 



n +( z *) = T3 d P* d Py / d P z 

fl J —oc J— CO J— oo 



exp-E&M* = z*~ —),Px,P y ,Pz]}- (29) 
/ m 



This condition yields 



)l -2|)-I J = 2, I . (30) 



The second condition expresses position independence of flow of momentum across the 
system. Equation ( p4[) is used with p z as the flow defining generating function. Condition 
dJ Pz /dz = yields 

* = §T ™ 

For flow of energy (heat) through the system, the generating function in equation (E3) is 
J2{p 2 /2m). If internal rotation is superimposed on translation (Eucken's correction [|33| ,|3?H ) , 
the relevant contribution should be added to the generating function. With atomic gases 
however, the result reads 

8 p z D pm 
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Heat conductivity (A) is the ratio between the sum of the exchanges at either walls (2 Jg) 
and the temperature gradient. Hence 

A = k B n (33) 
4 (3m 

In equations ([27|) and ([£]), the transport coefficients are expressed in terms of the ef- 
fective collision periodicity r. For absolute comparison with the experiment, an additional 
expression is required that relates the latter to the mechanical properties of the colliding 
species (mass and cross-section) at the given temperature. Without this additional in- 
formation, only the ratio between viscosity and heat conductivity may be compared with 
experimental data. This ratio is by definition Prandtl's number 

Pr = ^, (34) 
m A 

with c p as the constant pressure heat capacity. Implementation with the results obtained 
above yields the experimental results identically. 

B. Mixtures of atomic gases 

Let the components of a given mixture be indexed A and B, where A points to the sub- 
stance with the higher mass. Each component may be considered as a separate subsystem, 
with its own thermodynamic variables, interacting simultaneously with the other one and 
with the environment. The generalized Massieu function being extensive, we have for the 
composite system 

M=M A + M B - (35) 
With dilute gases or gases interacting as hard spheres, the individual generalized Massieu 



functions are defined as in equation fl22). For each component separately the generating func- 
tions to be used are the same as for single-component gases (see tables [I] and [IT]), excepting 
for the requirement of indexing the relevant masses in the appropriate generating functions. 
In stationary or quasi-stationary conditions (see below), for exchangeable properties where 
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equilibrium between the subsystems prevails, the intensities are the same. In the examples 
treated below, that is the case for the temperature and its moments and for the intensities 
conjugate to collective motion. Intensities conjugate to not exchangeable properties and 
their respective moments will be indexed according to the particular component they refer 
to. 

The collision periodicity has been shown above to be an essential ingredient in the dy- 
namics. In multi-component systems, there is an average collision periodicity for each of the 
constituents (ta, Tb). It measures for each how long the relevant atoms move freely before 
being halted by the matrix formed by the remaining particles. 

In multi-component systems there are homogeneous and heterogeneous collisions. Their 
frequencies add up. The efficiency for exchange of momentum from a colliding atom to the 
local thermodynamic bath depends on the masses of the collision partners. When a heavy 
particle hits a light constituent of the thermodynamic bath, its path is less disturbed and 
less momentum is transferred than in the opposite case. 

We assume a particle with mass m x and linear momentum P hitting a stationary matrix 
particle with mass m 2 . If the exit path of the matrix particle forms an angle ip with the 
incident one, momentum transferred to the matrix equals 2\P\ cos(i/j)rri2/ '{mi + 7712). Hence, 
the relative transfer efficiency of heterogeneous collisions is Im^jijni + m^). For the total 
effective collision frequency of atoms of one sort with respect to the matrix (1/r), the latter 
coefficient is the appropriate scaling factor relating efficiency of heterogeneous collisions to 
homogeneous ones. 

With hard spheres, the traditional expression for the collision periodicity between iden- 



tical particles with collisional cross-section d is known to be [33 



5 ^71171(3 

T = 16 I^P- (36) 

Let now the collision cross-sections be respectively g^a, c?bb, <^4B- Using the efficiency pa- 
rameter defined above and adding for either constituents the homogeneous and the hetero- 
geneous contributions, the total effective collision periodicities read 
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5 [p (n A d\ A 2m B 2 lm A + m B \ 1 
t a = — \ - + : n B d AB J— , 37 

5 //? / 2m A 2 lm A + m B n B d BB \ 1 
Tfi = 77^ V ~ , n Ad AB \ -z + , • 38 

From here on it is advisable to replace n A by xn and % by (1 — x)n. 

With hard spheres we have d AB = (d AA + d BB )/2. It appears that the published data 
on the viscosity of mixtures of atomic gases is accurate enough to allow the heterogeneous 
hard sphere cross-section to be corrected by a factor e close to 1. 

1. Viscosity 

With Couette flow conditions, the intensities under direct control of the surroundings 
are the intensities conjugate to the particle numbers of either substances (a a, &b), the 
temperature (or better 0) and the intensity conjugate to the linear moment of shear velocity 
<j y (see table |). Three intensities need still to be determined, namely the quadratic moment 
of the temperature (or better 72) and the quadratic moments of the particle distributions 
for A and B (0 M , 9 2 , B )- 

The three additional relations required for completing the thermodynamic description of 
the system are of the same vein as those used for Couette flow in single component gases. 
For symmetry reasons, it is easy to show that flow of shear momentum is independent of the 
particular values of three missing intensities. The principles involved in their determination 
will be therefore skipped. 

Flow of momentum is supported by either components. For each, the contribution is 
given according to equation (0), where the generating function to be implemented as <j)j is 
p y . Integration yields 

J Pv = \- XTA + ( X - X ^ ( 39 ) 

The viscosity of the mixture is therefore 

Vmix = ~a \ XT A + (1 - x)t b ], (40) 
29 



where n/(3 is the total pressure. 

The result of equation ( PH) is plotted in figure [I] for a mixture of Xe in He. Experimental 
results at 291 K [|3l| are indicated on the same graph (experimental uncertainties ~ ±1%). 
Correction factor e for heterogeneous collisions may be estimated by fitting the curve to the 
experimental results. That obtained without the correction factor is displayed as a dotted 
curve. The fit performed on the ten mixtures of atomic gases leads to the values of e in 
the range 1.03 - 0.98, the highest values being for mixtures with a light component (He). 
Results for the ten mixtures of atomic gases are listed in |39| . 



FIG. 1. Predicted and experimental viscosity of a mixture of Xe in He at 271 K. The smooth 
curve is for e = 0.98, the dotted curve for e = 1. 



2. Diffusion and thermal conduction 

For a binary system enclosed between two reservoirs at different temperatures separated 
by a distance 2D (x: mole fraction of the heavier substance (^4)), stationary conditions are 
completely described by the set of constraints listed in table|I|, the intensities being indexed 
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accordingly. 

The intensities conjugate to the particle numbers of either substances (a a, ozb), the 
temperature (or better (3) and its gradient (or better 71) are under direct control of the sur- 
roundings. Thermal interaction between the subsystems removes the necessity for indexing 
the latter two intensities. 

Three intensities remain to be determined, namely the two gradients of the particle 
distributions (9i t A, 9i,b) and the intensity conjugate to collective momentum from wall to 
wall (er z ). Hence, three additional conditions or equations are required. Two are identical to 
the conditions discussed for single component systems: mechanical equilibrium and global 
stationarity. 

Mechanical equilibrium of the system between its walls implies vanishing total pressure 
gradient. It does not require per se vanishing partial pressure gradient for either substances 
separately. A possible pressure gradient of A is indeed neutralized by an opposite gradient 
for B. By stating that the sum of the contributions of either substances to flow of momentum 
between the boundaries is position independent, the following equation is derived (compare 
equation fl3"Ip) 



The condition for global stationarity is defined along the same lines as above (equa- 
tions (p8|)— fl3~0f) )) where the density n + (z*) is now understood as the sum for the two compo- 
nents. As a result, the relation for internal collective motion (a z ) becomes (see equation (|30|)) 



The last condition to be considered concerns mutual diffusion or motion of the subsystems 



<f>j = 1, the particle flow of either subsystems is obtained. According to whether the named 
parameters are indexed A or B, the results are 




(41) 




(42) 



with respect to each other. By implementing equation (|24| ) with the generating function 



31 



The first contribution in either equations (that proportional to a z ) represents collective 
drag generated in the fluid by correlated effect of the walls. This acts on the two subsystems 
alike. Therefore it does not drive diffusion of the subsystems with respect to each other. By 
contrast, diffusion is related to the second part of the flow equations. As it may be verified, 
this is driven by the relevant partial pressure gradients. 

Diffusive stationarity is reached by differential displacement of the subsystems with re- 
spect to each other. Then we have for either subsystems vanishing partial pressure gradients. 
This is the remaining constraint for complete thermodynamic description of the system. In 
experimental conditions, the question is however whether diffusive stationarity has been 
reached in practical cases when thermal conductivity of multi-component mixtures is mea- 
sured. 

For measuring thermal conductivity, an appropriate binary mixture is prepared in a con- 
ventional thermostat. The walls are then brought at different temperatures. When possible 
acoustic perturbations have relaxed, the total pressure distribution is flat (equation (pEip). 
Nevertheless, on establishing the temperature gradient, pressure gradients of the individual 
constituents are created, forcing the particles to segregate. If the mixture consists of parti- 
cles with different mobility, final stationary conditions are slow to reach. The slower moving 
particles tend to remain distributed homogeneously, as they were before the temperature 
gradient was created. The partial pressure gradient of the faster moving subsystem com- 
pensates for resulting unbalance. Final equilibrium requires that the migration of the slower 
particles have taken place, cancelling all partial pressure gradients. 

The published experimental data that have been considered do not mention whether 
(or how much) the system has been allowed to relax the initially created individual partial 
pressure gradients. Let us assume this would not have occurred at all. The state is then 
pseudo-stationary as it continues to change slowly in time while the constituents still migrate 
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with respect to each other. The two subsystems should then be considered as acting inde- 
pendently for all the properties concerning the particle distributions. They remain however 
tightly coupled for the properties that are promptly interchanged. In particular, they share 
the same (3 and 71. The intensity a z conjugate to the collective momentum generated by 
the temperature gradient is also common to the two subsystems. Concerning the latter, its 
relation to the other intensities and to the collision periodicities is given by (|30D . Instead of 
( fl2f) we have now two relations, namely 

{oiA-n) ^to= 2 ^ (44) 
Kb - n ) ^ = 

By combining ( f4T|) with the two latter ones, an expression for the gradients of the indi- 
vidual partial pressures may be derived. Writing 

fl=^, (45) 

T B /m B 

this relation reads 

6l ' A ~2j-2 (l-x)R + x J' (46) 

In practical cases, depending on the mixture to be considered, when heat conductivity 
is measured, the system may be somewhere between the two extreme conditions. The 
uncertainty concerning how close diffusion has reached stationarity in the experimental 
conditions where the measurements have been performed, explains why thermal conduction 
data of mixtures are difficult to reproduce. Let us express the uncertainty by a coefficient 
c, to multiply the right-hand side of equation (f46|) . When discussing a homogeneous set 
of data with varying compositions x, we assume for simplicity that the same coefficient is 
valid. Equilibrium for diffusion implies c = 0. 
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FIG. 2. Thermal conductivity of a mixture of Xe in He at 271 K. The smooth curve is for c = 0.5 



Transport of heat is supported by either components of the mixture. For each, the con- 
tribution is given according to equation (|32"D, where the relevant intensities are determined 
as above. Hence, 
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(47) 



In comparing the result with experimental data, coefficient c may be taken as an ad- 
justable parameter. Using the data published by E. Thornton and coworkers jE]] for 291 K, 
coefficient c has been found to range between for light-light mixtures (e.g. He-Ne) and 1 
for heavy-heavy mixtures (Kr-Xe). Figure is an illustration of the results for the He-Xe 
mixture, where the value of c optimizes at 0.5. Results for the ten mixtures of atomic gases 
are listed in [55f|. Accuracy is better than the announced experimental precision of 4%. 



VI. STRUCTURE FORMATION 

A. Benard-Rayleigh thermal convection 

Stable vortices developing in a fluid bound by two horizontal plates at different tempera- 
tures in a vertical (gravitational) force field are named after Bnard and Rayleigh. Since their 
first description in 1900 they have been the subject of an abundant literature, being a typical 
example of structure formation in dynamic systems. The present discussion aims at devel- 
oping a set of differential equations based on Liouville's fundamental equation (^), where 
Hamiltonian iso-entropic contributions and stochastic irreversible interactions are clearly 
separated. 

When convection develops in a fluid, the simple vertical symmetry defined by the two 
temperature reservoirs and the external field acting perpendicularly to the walls is broken. 
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From unidimensional the problem becomes bi- or tridimensional, with increased mathemat- 
ical complexity. Convective motion stimulated by thermal strain may show many different 
patterns: rolls, cells, etc. The particular shape adopted by the system depends primarily on 
lateral boundary conditions. Depending on symmetry, we may want to have the coordinate 
system transformed. For simplicity we consider generation of rolls, for which Cartesian co- 
ordinates are most appropriate. Intensities referring to the different directions in space will 
be indexed accordingly. 

Let the fluid be confined between two parallel plates (distance 2D, taken as the 
z— direction) at different temperatures. The plates represent a double temperature reservoir 
defining a given average inverse temperature (3 and a gradient j z . 

The temperature gradient (7^) generates a linear moment of the density distribution. 
This causes a non- vanishing value of 9 Z . The relation between 9 Z and 7 2 depends on con- 
servation of momentum (compare (|3~ID). In the presence of gravity, pressure gradient is 



balanced by the external force (—mg), whence [p~5|j: 



0,-~j = -PmgD. (48) 

In field-free conditions, if 7 Z = 0, the system's center of mass is located half-way between 
the plates (Z = 0). Thermal strain and the external force displace the center of mass 
with respect to this neutral position. Expressed as a function of the intensive variables and 
using (f48|), the vertical moment of the particle density (N z ) is given by (|18l): 



N, = j(jt-/3mgD). (49) 

Gravitation orients the gradient of the particle density to the bottom regions but, if the 
system is heated from below, the distribution may reverse. This situation presents much 
analogy with the population inversion occurring in laser physics and the conclusions devel- 
oped below may readily be transferred to the domain of quantum optics [|4T~| , f4"2"|| . 

With convection, intricate distributions develop, requiring for their description additional 
or modified generating functions, supporting new intensive variables. The distributions will 
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be analyzed on the basis of generating functions analogous to the set listed in table |J, but 
extended and generalized to include the additional constraints. 

The geometric structure we want to focus on is periodic in the horizontal direction. The 
macroscopic wavelength is A, the rolls turning in alternate directions. Let the y— axis con- 
nect successive vortices, the x— direction being parallel to the motion's local global angular 
momentum. We isolate along the ^—direction a distance A/2 in which one vortex fits. This 
cell will be our system. 

Let us focus on the intensities conjugate to the vertical and horizontal moments of 
the particle density and investigate the role of the source/sink contribution in Liouville's 
equation ([231). When a vortex is active it perturbs the stationary distributions. The vertical 
gradients are modified and horizontal gradients develop. Let the vertical gradients dictated 
by the plates and valid in the absence of a vortex be indexed z, R from now on (e.g. Z ^). 
Because of slow relaxation of distortions due to finite transport in the fluid, in the presence 
of a vortex this changes to 9 Z . Thermal diffusivity, a process related to thermal conduction, 
tends to neutralize the latter change. Relaxation of horizontal gradients (indexed y) follows 
the same mechanism. The above may be summarized in the form of a set of linear dissipation 
equations: 

Z = —k z (9 z — 8 z r), 



(50) 

By hy Oy . 

with, considering the cell's geometry: 

(51) 



Hz — ( f) ) K ; 



.(*)+(*)■ 

(k: coefficient of thermal diffusivity). 

In a second step we investigate the conservative iso-entropic part of (|23|). 

Let N z and N y be the vertical and horizontal moments of the particle density. Using ( fj9|) 
and ( |48| ) it is clear that diV 2 = |d0 z . The same argument holds in the horizontal direction. 

By referring to equation flT5|), considering possible changes of the moments of the particle 



density, invariance of the entropy leads to: 
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dS = [3mgDdN z - 9 z dN z - 9 y dN y = 0. (52) 

where mgDdN z represents change of the potential energy as the system's centre of gravity 
moves vertically By combining the two first r.h.s. terms as 9* z = 9 Z — /3mgD, the constant 
entropy condition may be rewritten (6 1 *) 2 + 9y = constant. 

With u representing the angular velocity characterizing the vortex, iso-entropic circula- 
tion is described by the following set of differential equations: 

" (53) 

9 y = u9* . 

By combining the latter set with the equations for dissipation (|50D, the following set is 
obtained: 



9* = -u9 y - k z (9* z ~ 0t, R ) , 
6y = u Q*z ~ h®y 



(54) 



In laser physics, this set is named after Bloch. 

Gravitation acting on the horizontal density gradient activates the vortex while friction 
inhibits the collective motion. Complete balance between the conflicting forces has been 



elaborated in ||15|| . Full development cannot be given here. It may however be drafted by: 

2g G{X)u 

"n — 

5D y D 



^ -—0y ' ( 55 ) 



where function G(X) depends on the form factor of the vortex, v is the viscosity of the fluid. 
By combining equations ([54] ) and (|55|), threshold inversion conditions for vortex formation 
and the relevant form factor at and beyond threshold are easily determined (pitchfork bifur- 
cation). The results are in agreement with the experiment | 15| . Mutatis mutandis, in laser 
physics, the corresponding set of equations describes readily, next to the threshold require- 
ments, implications for bifurcations to a variety of unstable and chaotic working conditions. 



1. Von Karman type turbulence 



In fluid dynamics it is known that the dimensionless number named after Reynolds 
governs both turbulence in pipes or channels and the von Karman vortex streets produced 
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in flows past airfoils. This suggests that analogous mechanisms may be at work in the latter 
different kinds of systems. In order to stress the similarity, label von Karman might be 
generalized by attaching it to turbulent flows in ducts as well. 



The flows listed above have all been treated successfully elsewhere [[L5[] in the same 
thermodynamic context. The detailed mathematics being however rather cumbersome, only 
the principles involved will be outlined next. 

We consider a fluid flowing in a channel with constant width 2D (Poiseuille flow). In 
laminar conditions the velocity profile through the channel is known to be quadratic. Con- 
trasting with the Couette flow discussed above, the table of generating functions required 
for the thermodynamic description of the system implies therefore a quadratic function for 
the collective momentum contribution. For the same reason, the kinetic energy distribution 
involves a fourth order generating function. 

We are interested in the non-laminar regime. Instead of the above, let us suppose there- 
fore that a stream of vortices (with angular velocity uj) fitting exactly within the space be- 
tween the walls flows down the channel. The purpose is to examine this particular regime's 
stability. 

Contrasting with the classical procedure where the walls are taken as immobile, inves- 
tigating the vortices implies the observer to be running with the fluid at the same average 
speed {v y ). The walls are therefore taken to be moving in the opposite direction. 

The general dynamical equations for a vortex have been derived above (0). In the 
present case there is however no external force field to justify inversion. Symmetry breaking 
results however from the asymmetric field of collective kinetic energy prevailing between 
the channel's walls, caused by the very existence of the vortex superimposed on local av- 
erage downstream translational motion. Accurate analysis of the relevant distribution of 
downstream collective kinetic energy indicates that it contains at least a linear contribution 
(proportional to z, the direction perpendicular to the walls). Its gradient is proportional to 
the product v y uj. 

Following Bernoulli's theorem, a gradient of collective kinetic energy in a fluid generates 
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a pressure or a density gradient perpendicularly to the flow. This is equivalent to a force 
acting on the fluid, here in the z— direction. This contains the product v v uj. It depends 
also on the form factor of the relevant vortices. The modified set of dynamic equations may 
therefore be written: 

6t = -u6 y -k z [6t-C{X)Dv y u], 

(56) 

Vy = - kydy. 

To the latter set an equation for Cj must be added. 

The walls are the loci for mechanical interaction of the surroundings on the system. 
Confining the fluid within its boundaries represents a force exerted by the environment. 
Its magnitude is the local hydrodynamic pressure. If the system is symmetric with respect 
to facing boundaries, the forces exerted by the latter are equal and oppositely directed. 
In the present case, asymmetric collective kinetic energy distribution with respect to the 
boundaries (linear velocity superimposed on vortex) ensures unequal coupling of the fluid 
with either walls. 

When a particle collides with one of the walls, depending on whether this is the one where 
the local average velocity gradient is higher or lower, the reinjection trajectory following 
collision is more or less reoriented. Hence, the force acting by the walls on the system is 
tangent to the flow at one wall and perpendicular at the other wall. The perpendicular 
pressure forces at either walls are not antagonistic. The resultant (pressure = n/[3) acts 
mechanically on the system. If the particle density presents a gradient parallel to the walls 
(9 y 7^ 0) the resulting moment of the forces activates rotation of the vortex. Hence: 

W =4^ fl »-^- W - (57) 
In stationary conditions, the dotted functions vanish. By eliminating 9 Z and 6 y in the 
resulting set of equations, the following equation for uj results: 

k C 

uj 2 - Re uj + k y k z = . (58) 

oLrfjTTl 

where Re = 2Dv y / v is Reynold's number. 
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Besides the trivial solution (u = 0), expression d58|) is the characteristic equation for 
vortex stability. It has two roots (limit point bifurcation) . Depending on the value of the 
Reynolds number they may be real or complex. Only real solutions justify stable vortices. 
The value of the critical Reynolds number separating conditions for stable and unstable 
vortices depends on their form factor. The lowest value must be retained. It has been 
calculated in |15| for flow in a channel, in a square section pipe and past an airfoil. Results 
are in agreement with the experiment. 

VII. CONCLUSIONS 

The present paper reviews the state of affairs for solving the irreversibility paradox. 
The conflict raises from apparent contradiction between reversibility of microscopic laws of 
motion and the irreversible behaviour of macroscopic systems. The analysis convincingly 
demonstrates that attempts based on the assumption that the relevant macroscopic systems 
are perfectly isolated, cannot justify the law of increase of entropy. Such systems are indeed 
necessarily conservative, a property that holds for the entropy too. Some of the attempts to- 
wards escaping this fundamental symmetry property have been discussed. Their arguments 
have been shown to be inconclusive. As a result one is forced us to assume that perfect 
isolation conditions are incompatible with real physical systems. 

The notion of quasi-isolation has been introduced, indicating the condition of systems 
that are allowed to exchange energy fluctuations with their environment. With closed sys- 
tems, the environment is at least for a part represented by their walls. 

The conceptual basis for the theoretical investigation of the dynamics of quasi-isolated 
systems is furnished by the experimental evidences obtained from the analysis of Joule's 
experiment. 
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